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We use 300 K reflectance data to investigate the normal-state electrodynamics of the high tem- 
perature superconductor Bi2Sr2CaCu20s+i over a wide range of doping levels. The data show that 
at this temperature the free carriers are coupled to a continuous spectrum of fluctuations. Assuming 
the Marginal Fermi Liquid (MFL) form as a first approximation for the fluctuation spectrum, the 
doping-dependent coupling constant X(p) can be estimated directly from the slope of the reflectance 
spectrum. We find that X(p) decreases smoothly with the hole doping level, from underdoped sam- 
ples with p = 0.103 (T e = 67 K) where X(p) = 0.93 to overdoped samples with p = 0.226, (T c = 60 K) 
where X(p) = 0.53. An analysis of the intercept and curvature of the reflectance spectrum shows 
deviations from the MFL spectrum symmetrically placed at the optimal doping point p = 0.16. The 
Kubo formula for the conductivity gives a better fit to the experiments with the MFL spectrum 
up to 2000 cm' 1 and with an additional Drude component or an additional Lorentz component 
up to 7000 cm" 1 . By comparing three different model fits we conclude that the MFL channel is 
necessary for a good fit to the reflectance data. Finally, we note that the monotonic variation of the 
reflectance slope with doping provides us with an independent measure of the doping level for the 
Bi-2212 system. 

PACS numbers: 74.25.Gz, 74.62. Dh, 74.72.Hs 



I. INTRODUCTION 

The complete phase diagram of the high tempera- 
ture superconducting (HTSC) cuprates is not yet known. 
While there is universal agreement on the presence of an 
antifcrromagnetic phase at low doping and a supercon- 
ducting region with a maximum T c at a doping level of 
p = 0.16 holes per copper site at higher doping, there is 
considerable uncertainty about the nature of the normal 
state outside the superconducting phase, particularly in 
the overdoped region. It has been clear from the be- 
ginning of the study of HTSC materials that their nor- 
mal state is highly anomalous and not at all like the 
Fermi-liquid state of the conventional superconductors. 
Three features stand out. First, at high temperature, 
the resistivity and the optical reflectance are dominated 
by scattering processes with a linear variation of the 
scattering rate with frequency and temperature where 
1/r ~ max(w, 7rT). 1-3 Experimentally, these processes 
manifest themselves as a dc resistivity that varies linearly 
with temperature with a zero intercept on the tempera- 
ture axis 4 and an infrared reflectance that varies linearly 
with frequency over a very wide range of frequencies, well 
into to the midinfrared. 5 ' 6 This singular behavior of the 
scattering rate has been termed Marginal Fermi Liquid 
(MFL) behavior from the original postulate of Varma et 
al. 3 where it was assumed that the quasiparticle weight 
vanished logarithmically at the Fermi surface, and the 
Green's function of the carriers was entirely incoherent. 

The second anomalous feature of the normal state is 
the presence of a pseudogap, a depression in the den- 



sity of states below a temperature T* that decreases 
as the doping level p is increased. 7 The pseudogap re- 
gion on the phase diagram is not well defined and the 
exact value of T* depends on the experimental probe 
used. Photoemission and tunneling measurements of 
Bi 2 Sr 2 GaCu 2 8+5 (Bi-2212) point to a T* > 300 K for 
underdoped samples, 8 and below 300 K for optimally- 
doped samples. 8-11 In optical spectroscopy the pseudo- 
gap is seen as a depression of interplane (c-axis) conduc- 
tivity, which in underdoped YBCO with T c = 63 K, T* 
occurs below 300 K (Ref. 12). The NMR shift also be- 
gins to drop below its high temperature value near 300 
K (Ref. 13). In Bi-2212 angle-resolved photoemission 
(ARPES) measurements show that the pseudogap devel- 
ops below 250 K in slightly underdoped Bi-2212 (Ref. 14). 
Probably most accurate and recent results will be a c-axis 
transport measurement of Bi-2212 15 . In the literature T* 
is around 300 K at p = 0.16 with rather large error bars in 
the underdoped phases. In all, it appears that a temper- 
ature of 300 K is above T* for optimally and overdoped 
states and below T* for underdoped states in the Bi-2212 
system. 

The third normal-state anomaly is the development of 
a sharp scattering resonance in the a6-plane transport 
properties at low temperature. It manifests itself as a 
sharp increase in the a6-plane scattering rate at 500 cm_i 
that takes place below 150 K in underdoped materials 
and near T c in optimally doped ones. 16:17 It is also seen 
as a sharp kink in the ARPES dispersion curves. 14 

The aim of this work is to investigate, in a system- 
atic way, the transport properties in Bi-2212 for a wide 
range of doping levels through accurate reflectance spec- 
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troscopy at 300 K. We restrict our work to the high- 
temperature strange-metal region where the spectra are 
dominated by the MFL fluctuations avoiding the low 
temperatures near T c and where these fluctuations arc 
gapped by the appearance of the resonant mode. By con- 
fining our measurements to high temperatures, we also 
steer away from the problem of the phase diagram in 
the overdoped state. We have chosen Bi-2212 for our 
measurements for several reasons. First, this material 
has been widely investigated by the angle resolved pho- 
toemission and tunneling spectroscopy. Secondly, large 
single crystals can be grown and annealed to yield sam- 
ples spanning the overdoped as well as the undcrdopcd 
regions of the phase diagram. 

Previous work on the optical properties of Bi-2212 has 
concentrated on the optimally doped materials. 16,18-22 
The focus has been on analyzing the optical conductiv- 
ities or the scattering rates. Kubo-MFL analysis (see 
Sec. H-D for a detail description) has been used by Lit- 
tlewood et al., 23 Baraduc et al., 20 and Abrahams. 24 To 
improve the overall accuracy, we fit the reflectance curves 
to models, avoiding the additional Kramers-Kronig (K- 
K) analysis required to calculate the conductivity and 
the scattering rates (see Sec. II-D). We also produce the 
optical conductivities by using K-K analysis fit them and 
compare the results with those of reflectance fits. Since 
the reflectance is related to the scattering rate in a di- 
rect way, models where the scattering rate has a simple 
form, such as the MFL hypothesis, allow us to estimate 
the coupling constant X(p) directly from the reflectance 
data. 

The present paper is organized as follows. First, we 
introduce a simple way to extract approximate MFL pa- 
rameters directly from the reflectance data. Secondly, 
we report indications of deviations from the simple MFL 
form from the study of accurate 300 K (Bi-2212) re- 
flectance data in a wide range of doping levels. The data 
presented here focus on the overdoped region where we 
present new data for several highly overdoped samples 
(T c = 82, 73, 65 and 60 K). We also included data from 
earlier measurements, 16 ' 25 as well as work from Tu et al. 22 
As mentioned above, we focus on 300 K reflectance data 
for our analysis, first because we want to avoid the region 
where the sharp scattering resonance distorts the spectra 
and secondly because with our measurement technique 
the data at room temperature are the most accurate. 



II. EXPERIMENTAL DATA AND ANALYSIS 

A. Experimental Techniques 

We used Fourier-transform infrared (FTIR) spec- 
troscopy to obtain the reflectance data of floating zone 
grown single twinned crystals of Bi-2212. To obtain over- 
doped samples the crystals were annealed in 3 kbar liquid 
oxygen in sealed containers. 26 A polished stainless steel 
mirror was used as an intermediate reference to correct 



for instrumental drifts with time and temperature. An 
in-situ evaporated gold film on the sample was the final 
reflectance reference. 27 The reflectance of the gold films 
was in turn calibrated with a polished stainless steel sam- 
ple where we relied on Drude theory and the dc resistivity 
as the ultimate reference. An advantage of this technique 
is that it corrects for geometrical effects of an irregular 
surface. The in-situ gold evaporation technique gives ac- 
curate, better than ±0.5 %, room temperature data. 



B. Trends in Reflectance 

Fig. 1 shows the overall trends of the absorption (A = 
1 — R) with doping. Three representative samples are 
shown: an underdoped (T c = 67 K) sample, an opti- 
mally doped (T c = 91 K) one, and finally an overdoped 
(T c = 65 K) sample. We note that the absorption varies 
linearly with frequency at high frequency for all doping 
levels and that both the slopes of the curves and their 
intercepts with the absorption axis decrease with increas- 
ing doping. We also note that the optimally doped re- 
flectance shows best linearity up to high frequency (sec 
Fig. 4) but still we can use this universal linear trend in 
the relaxation region as a first approximation. The actual 
reflectance should have an additional second order doping 
dependent term, which we discuss with curvature analy- 
sis of the reflectance in the end of this subsection. Within 
the Drude model it is easy to show that in the relaxation 
regime (1/t <C ui <C u> p ), where the imaginary part of 
the refractive index is much larger than the real part, 
the absorbance A is given by A — 2/uj p t (Ref. 28), where 
1/t is the scattering rate, u> the frequency, and lo p the 
plasma frequency. Thus the linear variation of absorption 
suggests a linear variation of scattering rate 1 /r with fre- 
quency. We also show least-squares fits to straight lines 
in the frequency range from 500 to 1750 cm -1 in the 
figure. We will call this analysis "1 — R slope analysis" . 

At low frequency, below the relaxation region where 
lu <C 1/t, the absorption drops below the fitted lines. 
This is the Hagen-Rubens region where the reflectance 
varies as ^Juj. We note that the Hagen-Rubens fre- 
quency range gets smaller as the doping level increases, 
which confirms the notion that the scattering rate 1/t 
decreases with doping, since within the Drude model the 
changeover from Hagen-Rubens region to the relaxation 
region occurs at ujt = 1. 

We determined the doping level p from T c using the 
parabolic expression of Prcsland et al. 29 The expression 
is p(T c ) = 0.16 =F [1/82.6 (1 - TJT^)] 1 / 2 , where T™ ax 
is the maximum T c or T c (p = 0.16) and =F means that we 
use — for undcrdopcd samples and + for overdoped sam- 
ples. The determination of T™ ax is a delicate problem 30 
and, in the absence of a better method, we use the gener- 
ally accepted value of 91 K as the T™ ax for Bi-2212. We 
should mention here that we have one optimally-doped 
sample which is doped with additional small amount of 
Y to yield a relatively well ordered system and shows 
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FIG. 1: Absorption (1 - R) data at 300 K for 
Bi2Sr2CaCu2 0s+« at three representative doping levels: one 
underdoped (T c = 67 K) sample, one optimally-doped (T c = 
91 K) sample, and one overdoped (T c = 65 K) sample and 
their linear fits, shown as solid lines. We note that both the 
slope and the intercept of 1 — R decrease as the doping level 
increases. 
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FIG. 2: Doping dependent slope of the absorption curves 
shown in Fig. 1, S(p), of BiaSraCaCuaOg+a a 300 K. We ob- 
serve two different trends: a rapid drop of slope with doping 
in the underdoped region and a slower linear variation in the 
overdoped region with a crossover point at p = p opt = 0.16. 
The T c vs. p curve shows where our systems are in the phase 
diagram. 



a surprising ly high T c = 96 K 30 . The disadvantage of 
the Presland method is that it does not uniquely deter- 
mine the doping level of the sample since there are two 
independent p values for each value of T c . However, as 
Fig. 1 shows, the reflectance slopes vary monotonically 
with doping and provide an alternate method of deter- 
mining the doping level that does not suffer from this 
ambiguity. 

As a first step, we fit the 1 — R data to a straight line in 
the frequency range between 500 and 1750 cm -1 to get 
the slopes and intercepts at various doping levels. Fig. 2 
shows the doping dependence of the slope, S(p) obtained 
from the 1 — R curves this way. We note the smooth 
variation of the slope as the doping level changes. Two 
different trends can be discerned: a rapid decrease in the 
underdoped region and a slower decreasing trend in the 
overdoped region. 

Fig. 3 shows a doping dependent intercept, I(p) of 1— R 
on the lo = axis. The intercept is related to the scatter- 
ing rate at lu = and to the dc resistivity. The intercept 
also changes smoothly as the doping level changes. The 
intercepts are consistent with the slopes within the MFL 
hypothesis, except at the high doping range, where we 
observe a hint of a crossover around p — 0.21. There are 



three regions: a sharp decaying trend in the underdoped 
region similar to what we observed in the slope, a slower 
decrease in the optimally doped region similar to what we 
observed in the variation of the slope. At p > 0.21 we see 
a crossover where is the intercept drops more abruptly, 
i.e., the scattering rate drops quickly. 

We found another interesting doping dependent trend: 
a slight deviation from a straight line reflectance in the 
relaxation region giving rise to a curvature of the re- 
flectance data. We analyze the curvatures of 1 — R at var- 
ious doping levels in the following way: We take two fre- 
quency segments, 500—1250 enr 1 and 1500 — 2250 cm -1 
in the relaxation region. Then we fit each segment to a 
straight line. Next we get slopes and intercepts of the 
two fitted straight lines. We label these Sl and II (of 
the lower frequency section) and Sh and Ih (of the higher 
frequency section). We estimate two curvatures indepen- 
dently for each 1 — R curve: Cs = {Sl — Sh)/Sl from 
the slopes and Cj = — (Jz — Ih)/Il from the intercepts. 
Fig. 4 shows the two doping dependent curvatures. A 
positive (negative) curvature stands for an upward con- 
vex (concave) curve. As the figure shows, the curvature 
of 1 — R is positive for the underdoped samples, negative 
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FIG. 3: Doping dependent intercept, I(p), of 
Bi 2 Sr 2 CaCu2 8+ a at 300 K from the fits (see Table II). 
We observe three different trends: an abrupt decay in the 
underdoped region, a slow drop between p — p opt = 0.16 
and p = p c = 0.21, and an abrupt drop above p = p c as the 
doping level increases. We observe two crossover points at 
P = Popt and p = p c . We also show the T c vs. p curve. 



for the overdoped samples and goes through zero exactly 
at optimal doping where 1 — R is a straight line. 



C. Analysis of 1 — R 

The trends which we observe in Fig. 1, the linear varia- 
tion of the scattering rate with frequency and the overall 
decrease in the slopes and intercepts with doping, can 
be related to the parameters of the MFL theory where 
the scattering rate, 1/t(uj,p) w \(p)[u> + 7rT], where A is 
the coupling constant, p is a doping level, and T is the 
temperature. 23 

For thick enough superconducting samples the ab- 
sorbance (A) is 1 — i?. In the relaxation regime we have 
the following: 



A(uj, p ) 



2X s {p) , 2X i ( P )ttT 



w p (p)t(w,p) uj p {p) 



■U) 



u p (p) 



(1) 



The above equation is linear in frequency for a fixed 
doping level. We get the doping dependent slope, S(p), 
and intercept, I(p), from a least square fit of our 1 — R 



FIG. 4: Two doping dependent curvature of 1 — R of 
Bi2Sr2CaCu208+f at 300 K. The curvatures of the 1 — R 
curves go through zero at optimal doping where p = 0.16 (see 
Sec. II-B for a detailed description). 



data to a straight line. According to Eq. (1) the slopes 
and intercepts are related to each other: ttTS(p) = I(p) 
and therefore we would obtain the same value of the 
doping dependent coupling constant X(p) from either the 
slopes or the intercepts. To account for deviations from 
Eq. (1) in our approximate analysis we allow the two 
constants to be different. We call the coupling con- 
stants from the slope Ag(p) and from the intercept A/(p). 
Within the MFL hypothesis Ag(p) — A/(p). In terms of 
the measured quantities S(p) and I(p) they are given by: 



^s(p) = -uJ p (p)S(p), 



A/(p) = 



I{p)uj p (p) 
2ttT 



(2) 



To get the absolute value of coupling constants ui p (p) 
must be known. We estimate u) p (p) from the spectral 
weight of the conductivity below the interband-transition 
region. To determine the frequency of the onset of in- 
tcrband transitions, we note that the absorption coeffi- 
cients, a(u>) = ATra(oj)/[n(uj)c], begin to deviate from the 
low frequency form at 9000 cm -1 and extrapolate to zero 
at the same frequency (13000 cm -1 ) for all the samples, 
independent of doping, as shown in the upper panel of 
Fig. 5. This fact suggests the following method of esti- 
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FIG. 5: (a) Absorption coefficient a at 300 K for 
Bi2Sr2CaCu208+f at six representative doping levels: two 
underdoped (T c = 67 and 82 K), one optimally doped (T c = 
96 K), and two overdoped (T c = 82, 65, and 60 K) samples, 
(b) Squares of the doping-dependent plasma frequencies ujp 
of Bi-2212 at 300 K, calculated by using Eq. (3). The plasma 
frequency increases monotonically with doping. 



mating the spectral weight up to the interband transition: 

,.9000 i /-13000 



120 

7T 



i) 2 J 9000 



(3) 



In the lower panel of Fig. 5 we show the resulting 
plasma frequencies of Bi-2212 at 300 K obtained this 
way. We see that the plasma frequency increases mono- 
tonically as the doping level increases. 31 

Fig. 6 shows the doping-dependent coupling constants, 
Xs(p) and A/(p) obtained using Eq. (2). Since the 
doping-dependent plasma frequency increases smoothly 
with the doping level we see the same trends in Xs(p) and 
Xi(p) that we observed in the doping-dependent slope 
and intercept. We note that the coupling constant A/(p) 
is larger than Xs(p). This is not surprising since in our 
simple analysis we have neglected the frequency depen- 
dence of the plasma frequency required by causality. We 
correct this in the more complete analysis in the next 
section. The simple analysis shows an overall trend of a 
decrease of the dimensionless coupling constants as the 
doping level increases, which is consistent with the angle- 
resolved photoemission results. 32 
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FIG. 6: Doping-dependent coupling constants, Xs(p) and 
A/(p), of Bi 2 Sr 2 CaCu2 8+5 at 300 K. 



TABLE I: Kubo-MFL fit parameters of five representative 
Bi2Sr2CaCu208+i samples: two underdoped (T c — 67 and 
82 K), one optimally doped (T c —96 K), and two overdoped 
(T c = 82 and 65 K) samples. The results are shown in Fig. 7. 
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D. Kubo formula in MFL and Kubo-MFL fit 

1. Reflectance and Kubo-MFL fit 

Here we introduce a more accurate method of analysis. 
One can derive the complex optical conductivity from the 
Kubo formula and the MFL hypothesis as follows. 23,24 
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FIG. 7: Kubo-MFL fits at 300 K for Bi 2 Sr 2 CaCu20 8+ i at five 
representative doping levels: two underdoped (T c — 67 and 
82 K), one optimally doped (T c = 96 K), and two overdoped 
(T c — 82 and 65 K) samples. Symbols are fits and solid 
lines are corresponding data. The fit parameters are shown 
in Table I. 



where 



'E R (ui) = Xui In i — x, 

lu c 2 

Z a (lu) = Z r (lu)*, 

x = max (| a; |, ttT). 



(5) 



Here T, R (co) and T, A (tu) are the retarded and advanced 
single-particle self-energies, respectively, ui c is the cutoff 
frequency and (3 = 1/kgT. 
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where e(u>) is the complex dielectric function, R(uj) is the 
reflectance and e w is the background dielectric constant. 
We will refer to a fit using Eqs. (4), (5) and (6) as a 
"Kubo-MFL fit". We fixed the plasma frequencies, 
and the cut-off frequency, leaving only A as a free param- 
eter for each doping level. 

Fig. 7 shows five representative data sets [T c — 67 K 
(UD), T c = 82 K (UD), T c = 96 K (OPT), T c = 82 K 
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FIG. 8: Kubo-MFL fits for extended frequency range at 300 K 
for Bi2Sr2CaCu208+i at five representative doping levels: 
two underdoped (T c = 70 and 82 K), one optimally doped 
(T c = 96 K), and two overdoped (T c — 82 and 65 K) samples. 
The symbols (MFL+Drude), dotted lines (MFL+Lorentz) 
and dash-doulbledotted lines (Drude+Lorentz) are the fits 
and the solid lines are the corresponding data. The fit pa- 
rameters are shown in Table III. The inset shows a magnified 
view of the fits and demonstrates the need for an MFL con- 
tribution for a good fit to the data. 



(OD), and T c = 65 K (OD)] and their Kubo-MFL fits 
between 60 and 1700 cm -1 and Table I shows param- 
eters for the fits. One interesting feature of the fits is 
that the fitted curves have slight curvature while the re- 
flectance data are straighter. The curvature of the fits 
above 500 cm -1 is convex upward (see Sec. II-B) sup- 
porting the argument above. The pseudogap induces a 
concave upward curvature on the reflectance. 

Table II shows the slope and intercept for each Bi-2212 
sample as well as the coupling constant from the 1 — R 
slope analysis between 500 and 1750 cm -1 for each dop- 
ing level, as well as the coupling constants from Kubo- 
MFL fits between 60 and 1700 cm -1 , shown in the Xk{p) 
column. 

We were unable to get a fit to the data up to 7000 cm -1 
with the simple MFL parameterization of the data. To 
get a reasonable fit for an extended range of frequencies 
we have to add a parallel Drude channel to the conduc- 
tivity. We fit the same Bi-2212 data at 300 K from 60 to 
7000 cm -1 by adding one Drude oscillator. Here we have 
three free parameters: the coupling constant A, width 
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FIG. 9: Kubo-MFL fits for extended frequency range at 300 
K for Bi2Sr2CaCu2C>8+<5 optical conductivity at five represen- 
tative doping levels: two underdoped (T c — 70 and 82 K), one 
optimally doped (T c = 96 K), and two overdoped (T c = 82 
and 65 K) samples. The symbols (MFL+Drude) and the dot- 
ted lines (MFL+Lorentz) are fits and the solid lines are the 
corresponding data. The fit parameters are shown in Ta- 
ble III. 



and strength of the Drude oscillator 70 and uj p d ■ Fig. 8 
shows the data and Kubo-MFL fits for the extended spec- 
tral range and Table III shows parameters for the fits. We 
observe that the coupling constants [Xk d {p)] are slightly 
higher than those [Aif (p)] of the previous Kubo-MFL fits 
without oscillators. The spectral weight of the Drude os- 
cillator needed for the good fit shown in the figure ranges 
from 44 % of the MFL spectral weight for the underdoped 
sample down to 17 % for the highly overdoped sample. 
The monotonic variation of the added Drude component 
with doping suggests that the deviations from the MFL 
form at high frequency may have a physical significance. 

We also fit the same data between 60 and 7000 cm -1 
with the MFL channel and one Lorentz oscillator. Here 
we have four free parameters: the coupling constant A, 
center frequency, width and strength of the Lorentz oscil- 
lator Ui, 7l and u> p l- Fig. 8 shows the data and Kubo- 
MFL fits for the extended spectral range and Table IV 
shows parameters for the fits. The spectral weight of the 
Lorentz oscillator needed for the good fit shown in the 
figure ranges from 38 % of the MFL spectral weight for 
the underdoped sample down to 13 % for the highly over- 
doped sample. We also calculate the mean square devia- 



FIG. 10: Comparison between two fits: MFL channel + a 
Drude channel and MFL channel + a Lorentz channel for 
reflectance and optical conductivity data. 



tion x 2 of the two high frequency fits to compare the qual- 
ities of the fits. The results are shown in Table VI. With 
one additional fit parameter the MFL+Lorentz yields a 
slight improvement in the fit. The exact physical mean- 
ing of the additional Drude or Lorentz oscillator is not 
clear but deviations from the simple MFL parametriza- 
tion is expected once the energies begin to approach the 
widths of the bands. In Fig. 10 we show the resulting cou- 
pling constants, Xkd(p) with Xkl{p) for both reflectance 
and the optical conductivity (see the following subsec- 
tion) . 

We fit the same reflectance data between 60 and 7000 
cm -1 with a Drude channel and one Lorentz channel. 
The fit parameters and resulting fits are shown in Ta- 
ble V and Fig. 8, respectively. Interestingly in the two 
overdoped samples (T c = 82 K and T c — 65 K) the least 
square fits converge to two Drude oscillators even though 
we start with one Drude and one Lorentz oscillator. We 
also show the mean square deviations Xdl m Table VI. 
Even though the Drude+Lorentz has the largest number 
of parameters it shows the worst fit. This suggests that 
an MFL channel is necessary for a good fit to the Bi-2212 
data and probably other cuprates as well. 
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TABLE II: The table shows the sources of 300 K 
Bi 2 Sr 2 CaCu2 08+5 reflectance (R) data. We fit the 1 — R 
data to a straight line in the frequency range between 500 
and 1750 cm" 1 . The slope and intercept of the fits are shown 
in the table. The sample whose T c is 96 K is in its optimal 
doping state. It was assigned a doping level of 0.16. We also 
show the coupling constants from Kubo-MFL fits between 60 
and 1700 cm -1 in the table. 



T C (K) 


V 


slope (/um) 


intercept 


As(p) 


Ax(p) 


Ref. 


67 


0.103 


1.069 


0.0818 


0.809 


0.899 


33 


70 


0.107 


1.232 


0.0884 


0.927 


1.165 


34 


82 


0.125 


0.844 


0.0670 


0.688 


0.720 


34 


89 


0.144 


0.796 


0.0620 


0.657 


0.658 




91 


0.160 


0.790 


0.0630 


0.670 


0.702 


22 


96 


0.160 


0.803 


0.0616 


0.681 


0.693 


new, 30 


90 


0.172 


0.669 


0.0616 


0.585 


0.606 


34 


82 


0.195 


0.625 


0.0555 


0.587 


0.606 


new,30 


73 


0.209 


0.594 


0.0550 


0.561 


0.584 


new 


70 


0.213 


0.552 


0.0520 


0.524 


0.533 


34 


65 


0.219 


0.537 


0.0418 


0.517 


0.457 


new 


60 


0.226 


0.564 


0.0470 


0.526 


0.506 


new 



TABLE IV: The MFL+Lorentz fit parameters of five repre- 
sentative Bi2Sr2CaCu20s+« samples: two underdoped (T c = 
67 and 82 K), one optimally doped (T c =96 K), and two over- 
doped {T c = 82 and 65 K) samples for both reflectance and 
conductivity. The results of the fits are shown in Fig. 8 and 
9. We have added to the MFL conductivity channel (plasma 
frequency u> p ) a Lorentz channel with a center frequency lul, 
a plasma frequency lu p l and a width 7l- In the table " (R)" 
and "(<ri)" stand for the reflectance and the conductivity, re- 
spectively. 



Tc (K) 


oj c (cm r ) 




Up (cm 1 ) 


too 




1L 


^>pL 


67 (R) 


7100 


0.884 


15040 


3.14 


1229 


5161 


9344 


82 (R) 


7100 


0.731 


16320 


3.21 


1003 


4600 


8910 


96 (R) 


7100 


0.706 


16980 


3.52 


801 


4177 


9614 


82 (R) 


7100 


0.612 


18789 


4.16 


741 


3470 


9851 


65 (R) 


7100 


0.461 


19270 


4.25 


1472 


3436 


6944 


67 (<7i ) 


7100 


0.977 


15040 


3.14 


373 


4115 


10765 


82 ((7i ) 


7100 


0.721 


16320 


3.21 


503 


3585 


10216 


96 (ai ) 


7100 


0.698 


16980 


3.52 


606 


3007 


10533 


82 (o-i) 


7100 


0.640 


18789 


4.16 


492 


2322 


10853 


65 (<7i ) 


7100 


0.447 


19270 


4.25 


463 


2722 


9555 



TABLE III: The MFL+Drude fit parameters of five represen- 
tative Bi2Sr2CaCu2 0s+« samples: two underdoped (T c = 67 
and 82 K), one optimally doped (T c =96 K), and two over- 
doped {T c — 82 and 65 K) samples for both reflectance and 
conductivity. The results of the fits are shown in Fig. 8 and 
9. We have added to the MFL conductivity channel (plasma 
frequency u p ) a Drude channel with a plasma frequency u) p d 
and a width jd- In the table "(7?)" and "(<7i)" stand for the 
reflectance and the conductivity, respectively. 



Tc (K) 


Lu c (cm r ) 


\kd 


u!p (cm 1 ) 


too 


7d 


(*JpD 


67 (R) 


7100 


1.019 


15040 


3.14 


4790 


9970 


82 (R) 


7100 


0.797 


16320 


3.21 


4250 


9360 


96 (R) 


7100 


0.749 


16980 


3.52 


4000 


9960 


82 (R) 


7100 


0.647 


18789 


4.16 


3300 


10260 


65 (R) 


7100 


0.506 


19270 


4.25 


2930 


7920 


67 (o-O 


7100 


1.001 


15040 


3.14 


3980 


10843 


82 (o-i) 


7100 


0.762 


16320 


3.21 


3270 


10397 


96 (<ti) 


7100 


0.789 


16980 


3.52 


2708 


10987 


82 (ai ) 


7100 


0.730 


18789 


4.16 


2013 


11400 


65 (ai ) 


7100 


0.461 


19270 


4.25 


2569 


9497 



2. Conductivity and Kubo-MFL fit 

We also fit the optical conductivity between 60 and 
7000 cm" 1 with an MFL channel and two different oscil- 
lators: MFL + Drude or MFL + Lorentz. We determined 
the optical conductivity from the measured reflectance by 
using a Kramers-Kronig analysis 35 , for which extrapola- 
tions to w — > and oo must be supplied. For lo — > 0, 
the reflectance was extrapolated by assuming a Hagen- 
Rubens frequency dependence, (1 — R) oc cj 1 / 2 . The re- 
flectance has been extended to high-frequency by using a 
literature data 36 and free-electron behavior (R ex w" 4 ). 
The result of the fits to the data are shown in Fig. 8 and 
the fit parameters are given in Table III and IV. We also 



TABLE V: The Drude+Lorentz fit parameters of five repre- 
sentative Bi2Sr2CaCu2C>8+i samples: two underdoped (T c = 
67 and 82 K), one optimally doped (T c =96 K), and two over- 
doped (T c = 82 and 65 K) samples for reflectance. The results 
of the fits are shown in Fig. 8. We have five fit parameters: 
two from a Drude channel with a plasma frequency oj p d and 
a width 7u and three from a Lorentz channel with a center 
frequency ujl, a plasma frequency lu v l and a width jl. 



T c (K) 


too 


7d 


u)pD (cm 1 ) 




JL 


<jJ p L 


67 (R) 


3.14 


698 


9149 


2336 


8762 


13200 


82 (R) 


3.21 


618 


10300 


2201 


8493 


13323 


96 (R) 


3.52 


607 


10810 


2022 


7343 


13187 


82 (R) 


4.16 


625 


11958 


0.14 


6950 


14735 


65 (R) 


4.25 


444 


12268 


0.12 


6040 


14491 



TABLE VI: Comparison between three fits for reflectance 
data: MFL channel + a Drude channel, MFL channel + 
a Lorentz channel and a Drude channel + a Lorentz chan- 
nel. We named the mean square deviations Xmd f° r t ne 
MFL+Drude fit, Xml for the MFL+Lorentz fit, and \dl 
for the Drude+Lorentz. The additional parameters of the 
Lorentz fit yields only a slight improvement when we com- 
pare the MFL+Lorentz fit with the MFL+Drude fit. 



Tc (K) 


Xmd 


Xml 


Xdl 


67 


18.65 


17.03 


21.45 


82 


16.35 


15.63 


18.62 


96 


2.21 


1.48 


6.59 


82 


3.61 


2.89 


6.14 


65 


7.99 


6.57 


43.70 


No. of parameters 


3 


4 


5 
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show the resulting doping dependent coupling constants 
from two different fits for both reflectance and conduc- 
tivity in Fig. 9. We observe smoother doping dependent 
coupling constant in the result from the reflectance fit. 
This is one of the reasons why we use reflectance to study 
the doping dependent properties. 

The additional Drude and Lorenz oscillators intro- 
duced to fit the wider frequency range data should not 
be taken literally as additional physical channels of con- 
ductivity. They only serve to correct the frequency de- 
pendence of the MFL spectrum that is known not to fit 
well over a wider range of frequencies 19,20 ' 37 . 



III. DISCUSSIONS 

We note that as a first approximation, the MFL spec- 
trum of fluctuations accounts for the data well in the 
low frequency region of u < 2000 cm -1 both within the 
simple 1 — R slope analysis and the more accurate Kubo 
formalism. In Fig. 11 we compare the dimensionless cou- 
pling constants X(p) from the two different methods of 
analysis: l-R and Kubo-MFL fit (without oscillators). 
The coupling constant decreases uniformly with doping 
except for the highly underdoped samples. However, in 
the low doping range the pseudogap temperature T* may 
approach our measurement temperature of 300 K. The 
pseudogap leads to a decreased scattering rate at low fre- 
quencies which would in turn lead to an enhanced slope 
of the absorption curve. The upturn in the slope at low 
doping levels may therefore well be a pseudogap effect. 
Thus, our first conclusion is that at room temperature the 
coupling constant X{p) decreases uniformly in the doping 
region that we have studied without any evidence of any 
crossovers. 

We also show in Fig. 11 the coupling constants from 
the 1 — R slope analysis and from ARPES. 32 Johnson et 
al. 32 argued that in the normal state the self-energy is 
well described by the MFL hypothesis and in the over- 
doped region the difference between the superconducting 
and normal state dispersion vanishes. They calculated 
the coupling constants by using A = — (<9Re£/i9w)£; F . 
We have augmented their analysis by estimating the cou- 
pling constant at room temperature which is lower than 
their published low-temperature value in the underdoped 
region. This is probably due to a rearrangement of spec- 
tral weight of the spin-fluctuation spectrum with temper- 
ature. As the temperature is lowered, this rearrangement 
takes two forms: spectral weight is removed from low fre- 
quencies as a result of the development of the spin gap 
but added to intermediate frequencies with the devel- 
opment of the neutron resonance that couples strongly 
to the carriers. The net result is a strong increase in A 
which suggests that spectral weight is removed from high 
frequencies to fill the neutron mode. 38 

The overall good agreement of the coupling constant 
and its doping dependence between the infrared and the 
ARPES data is surprising since the infrared response rep- 
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FIG. 11: Coupling constants, X(p), of Bi2Sr2CaCu20s+a at 
300 K from two different methods, discussed in the text and 
one result from angle resolved photoemission spectroscopy of 
Bi-2212 data in the superconducting state. 32 



resents an average over the Fermi surface, whereas the 
ARPES data shown in Fig. 11 represents corrections to 
the self energy in the (it, 7r) direction. Also, ARPES mea- 
sures the particle life-times directly, whereas infrared is 
weighted in favor of large angle scattering. 

Our data on the intercept and the curvature does show 
evidence of doping-dependent discontinuities. The effect 
of the pseudogap on the intercept would be opposite from 
that on the slope, causing it to be smaller at low doping 
if there was a pseudogap. Instead, as shown in Fig. 3, we 
observe the opposite effect in the underdoped region, the 
intercept is higher than what would be expected from a 
uniformly increasing A with underdoping. Also, in the 
highly overdoped region the intercept drops, whereas the 
slope appears to have a more uniform trend in the 0.16 to 
0.22 doping region. The overall result of these effects is to 
produce an S-shaped curve centered at optimal doping. 

However, the clearest evidence of a singular doping 
level at p = 0.16 comes from our analysis of the dop- 
ing dependence of the curvature plotted in Fig. 4. Here 
we see a good fit to a straight line that goes through zero 
exactly at optimal doping p = 0.16. This phenomenon 
is directly related to the well known observation that 
the dc resistivity in high temperature superconductors 
varies linearly with temperature only at optimal dop- 
ing and curved on either side of optimal doping. This 
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was shown by Takagi et al. 39 for overdoped and under- 
doped L^-zSr^CuO^ by Carrington et al. in under- 
doped YBa2Cu307_5 (Ref. 40), by Kubo et al. in over- 
doped Tl 2 Ba 2 Cu02 (Ref. 41) and by Konstantinovic et 
al. for undcrdoped Bi-22 1 2 42 and for underdoped and 
overdoped Bi2Sr1.6Lao.4CuO,, 43 . The present work ex- 
tends this work to Bi-2212 and covers both overdoped 
and underdoped materials for the same system. Needless 
to say, we study the frequency dependence whereas the 
dc resistivity work is based on the temperature depen- 
dence. Thus we show by an independent method that at 
room temperature the MFL fluctuations are centered at 
optimal doping p = 0.16. 

IV. CONCLUSIONS 

We have analyzed the 300 K reflectance data of Bi-2212 
single crystals at various doping levels from underdoped 
(p = 0.103) to highly overdoped (p = 0.226). We found 
three smoothly varying quantities with doping, the dop- 
ing dependent slope, intercept and the curvature of the 



reflectance in the relaxation region of frequencies. From 
the smooth trend of the reflectance slope with doping 
we derive a useful measure of the doping level. We es- 
timate a reliable doping dependent dimensionless cou- 
pling constant As (p) from the 1 — R slope analysis which 
varies smoothly from 0.93 to 0.52 in the doping range 
from p = 0.103 to 0.226. We also obtained doping de- 
pendent curvature of 1 — R. Plotting the curvature as 
a function of doping we found that the curvature goes 
through zero at optimal doping. 
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